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Young accreted globular clusters in the outer halo of M31 
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ABSTRACT 

We report on Gemini/GMOS observations of two newly discovered globular clusters in 
the outskirts of M31. These objects, PAndAS-7 and PAndAS-8, lie at a galactocentric 
radius of « 87 kpc and are projected, with separation ps 19 kpc, onto a field halo 
substructure known as the South- West Cloud. We measure radial velocities for the two 
clusters which confirm that they are almost certainly physically associated with this 
feature. Colour-magnitude diagrams reveal strikingly short, exclusively red horizontal 
branches in both PA-7 and PA-8; both also have photometric [Fe/H] = — 1.35±0.15. At 
this metallicity, the morphology of the horizontal branch is maximally sensitive to age, 
and we use the distinctive configurations seen in PA-7 and PA-8 to demonstrate that 
both objects are very likely to be at least 2 Gyr younger than the oldest Milky Way 
globular clusters. Our observations provide strong evidence for young globular clusters 
being accreted into the remote outer regions of M31 in a manner entirely consistent 
with the established picture for the Milky Way, and add credence to the idea that 
similar processes play a central role in determining the composition of globular cluster 
systems in large spiral galaxies in general. 

Key words: globular clusters: general, galaxies: individual: M31 
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1 INTRODUCTION 

It has long been recognized that some fraction of the glob- 
ular clusters in the Milky Way halo have very likely been 
accumulated via the accretion and destruction of their satel- 
lite host galaxies, rather than having formed in situ. How- 
ever much of the observational support for this picture is 
circumstantial rathe r than direct. In their seminal paper 
ISearle fc Zinnl |l978) showed that halo clusters outside the 
solar circle exhibit no abundance gradient with Galactocen- 
tric radius, and further that the horizontal branch (HB) 
morphologies for these objects (taken as a proxy for age) 
show little correlation with abundance. They postulated 
that the outer globular clusters did not form during the ini- 
tial collapse of the Galaxy but originated in smaller "proto- 
galactic fragments" that were subsequently accreted into the 
Galactic potential well, even long after the formation of the 
central regions of the Milky Way was complete. 

Modern studies with vastly improved datasets still 
support this scenario. It is now known that a wide 



variety of properties commonly associated with glob- 
ular clusters in the outer parts of the Milky Way, 
such as their ages (measured directly from precise 
main-sequence turn-off photometry), kinematics, HB mor- 
phologies, spatial locations, luminosities, and sizes , are 
consistent with an external origin (see e.g., Zinnl 
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that age is the dominant s econd param eter controlling hori- 
zontal branch morphology, I Zinnl (|l99ot ) labelled the ensem- 
bles of globular clusters with red and blue HBs at given 
[Fe/H] as the "young" and "old" halo systems, respectively. 
That young halo clusters typically have ages several Gyr 
younger than old halo members has now been explicitly 
demonstrated fe.g.. |Potter et al]|2010l . l201ll h and it is these 
objects, with their characteristic diffuse structures and fre- 
quently retrograde orbits, that are the strongest candidates 
for having been accreted into the Milky Way. 
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The "smoking gun" for the accretion scenario 
was the discovery of t he d isrupting Sagittarius dwarf 
lllbata. Gilmore fc Irwinl 1 19941 ) and the recognition that 
at least four remote Galactic globular clusters (M54, 
Arp 2, Terzan 7 and 8) are in fact members of 
this galaxy caught in t he midst of their arrival 
into the Milky Way halo Jlbata, Gilmore fc Irwinl 1 19951 ; 
Da Costa fc Armandrofj Il995h . It has been shown more 
recently that at least one young halo globular cluster 
widely separated from the main body of Sagittarius on 
the sky, Palomar 12, also on c e belonged t o that galaxy 
<|Martmez-Delgado et all |2002| ; ICohenl [2004); furthermore, 
a handful of additional clusters are strong candidates for 
being ex- members (e.g.. iBellazzini. Ferraro fc Ibatal 120031 ; 
lLaw fc Maiewskill2010l ; iForbes fc Bridgesll2010r i. 

If a significant fraction of remote Milky Way globular 
clusters have indeed been accreted, it might be expected that 
the outskirts of these objects would be an excellent place to 
search for the remnants of their host gala xies. However, de- 
spite a significant investment o f effort (e.g.. lSohn et al 1l2003l ; 
iMartmez-Delgado et alj 2004 ) and a few tantalising re cent 
hints (e.g.. lOlszewski et all 120091 ; ISollima et ail |20T2j ). no 
such streams have yet been definitively identified. To date, 
the handful of clusters arriving along with the Sagittarius 
dwarf represent the only direct evidence for the accretion 
scenario; all else remains circumstantial. In particular, de- 
spite the compelling properties of the young halo ensemble, 
it is not clear to what extent these objects are representa- 
tive of the entirety of the accreted population - several of 
the Sagittarius clusters are old halo members, and indeed a 
number of the clusters seen in the LMC a nd the Fornax dSph 
woul d also be classified as such (e. g., iMackev fc Gilmord 
|2004| ; iMackev fc van den Berghl^OoH ). 

M31 is of vital importance in providing a unique ex- 
ternal system for addressing the above problems; it is the 
only other large galaxy in which globular clusters may be 
resolved into individual stars. Wide-field surveys conducted 
by our group, including, most re cently, the Pan-Andromed a 
Archaeological Survey (PAndAS; lMcConnachie et al.| [2009). 
have uncovered a wealth of streams and substr u cture in the 
M31 field halo (see also iFerguson et al.l |2002| ; llbata et al] 
120071 ). along with a large number of remote globular clus- 
ters with projected radi i in the range m 15 — 145 kpc 
jHuxor et al.l 120051 , 120081 . 2013, in prep.). We have previ- 
ously demonstrated that there is a striking correlation be- 
tween a number of the tidal streams seen in the M31 halo and 
the po sitions of many of the globular clusters |Mackev et al.l 
l2010al ); now that the PAndAS footprint is complete, it has 
become clear that this association apparently holds across 
the entire M31 halo beyond w 30 kpc (Mackey et al. 2013, 
in prep.). The unavoidable conclusion from these observa- 
tions is that much of the outer M31 globular cluster system 
(> 80%) has been assembled via the accretion of cluster- 
bearing satellite galaxies, just as has long been suspected in 
the Milky Way. 

We have embarked on a detailed observational pro- 
gramme aimed at measuring (i) radial velocities for, and 
(ii) the resolved properties of, remote M31 globular clusters, 
especially those projected onto tidal streams and overdensi- 
ties. These observations will confirm (or refute) individual 
associations between clusters and substructures, and allow 
us to explore in detail the characteristics - in particular the 



metallicities, HB morphologies, and structures - of accreted 
M31 globular clusters for direct comparison with Milky Way 
subsystems such as the young and old halo ensembles. 

In this paper we report on some of the first results 
of this work. We have targeted two newly-discovered M31 
globular clusters from the PAndAS survey, PAndAS-7 and 
PAndAS-8 (Huxor et al. 2013, in prep.), with imaging and 
spectroscopic observations using the GMOS instrument on 
Gemini North. Assuming the usual M31 distance modulus of 
(m - M) = 24.47, PA-7 and PA-8 lie at projected galacto- 
centric radii of 86 and 88 kpc, respectively, and are separated 
by a distance of 81.7' ~ 18.6 kpc. Both project onto a field 
halo substructure kn own as the South- We st Cloud, as shown 
in Figure [1] (see also IMackev et al.|[2oToal '). 

So far little is known about either the target clusters 
or the South-West Cloud. Based on the PAndAS discov- 
ery images both PA-7 and PA-8 appear sub-luminous, with 
M v « -4.9 and -5.3, respectively. With (V - I) ~ 1.0, 
they are also mildly redder than other globular clusters at 
simila r radii (which t ypically have (V — I)o ~ 0.9 as in 
e.g jHuxor et alj|201l]). The South-W est Cloud was first de- 
scribed bv lMcConnachie et al l |2009), and is one of the more 
diffuse substructures identified in the M31 halo to date. Den- 
sity maps construc ted using different photo metric cuts (such 
as Fig. [T] see also I Richardson et al .1120111 ). indicate that it 
must be predominantly composed of stars with [Fe/H] < — 1. 
No kinematics h ave yet been measu red for the South-West 
Cloud; however, iLewis et al.l l|2012l 'l describe a large spur 
of Hi gas extending from the central parts of M31 to very 
nearly overlap with this stellar substructure. While any link 
between the two remains circumstantial, we note that the 
Hi velocity is ~ -470 kms -1 , with width « 200 kms -1 . 

In what follows we describe our data acquisition and re- 
duction (Section [2]), the results inferred from radial velocity 
measurements and the cluster colour-magnitude diagrams 
(Section [3]), and the implications for our understanding of 
accretion processes in assembling both M31 and the Milky 
Way globular cluster systems (Section gj. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Imaging and photometry 

We imaged PA-7 and PA-8 over a number of nights in Au- 
gust and September 2008, as indicated in Table [1] using 
the Gemini Multi-Object Spectrograph (GMOS) at the 8.1m 
Gemini North telescope on Mauna Kea, Hawaii. The data 
were obtained in queue mode via program GN-2008B-Q-22 
(PI: Mackey) during clear, photometric conditions and un- 
der excellent seeing (0.4" - 0.65"). 

The GMOS imager |Hook et alj|2004 ) is comprised of 
three adjacent 2048 x 4096 pixel CCDs separated by gaps 
of ~ 2.8", and has a field of view (which does not cover 
the full CCD package) of 5.5' x 5.5'. To take advantage of 
the high quality conditions we employed unbinned imaging, 
resulting in a plate scale of 0.0727 arcsec/pixel. We obtained 
our observations in the GMOS g 1 and i! filters, which are 
similar, but not identical, to the g and i filters used by the 
Sloan Digital Sky Survey (SDSS), and the g' and i! filters 
employed by the MegaCam imager on the Canada-France- 
Hawaii Telescope (CFHT). We took six images with the g' 
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Figure 1. Spatial density map, excised from the final PAndAS footprint, of stellar sources possessing luminosities and colours consistent 
with being metal-poor red giant branch stars ([Fe/H] < —1.4) in the M31 halo. North is up and east is to the left. The two dashed 
circles represent projected galactocentric radii of 30 and 130 kpc (2.2 and 9.6 degre es). Globular c lusters are marked with red points. 
In the inner regions of M31 these come from the Revised Bologna Catalogue llGalleti et alj|2004h ; outside £S 30 kpc almost all have 
been discovered from our wide-field surveys, including PAndAS. The two clusters on which we focus in this paper, PA-7 and PA-8, 
are labelled. These project onto the South- West Cloud, indicated with arrows. A third cluster, not studied here, also sits on this field 
substructure - this is PA-14, to the south-east of PA-8. Also not able in the field of view are the narrow Nort h- West Stream, the Giant 
Stream to the south, and two tangential streams to the east (see llbata et al. 2007; McCoimachic ct al. 20091). 



Table 1. Log of imaging and spectroscopic observations of the two target clusters and two velocity reference clusters. 



Cluster 


Coordinates (J2000.0) 


Filter or 


Number of 


Total 


Dates 


Image 


Name 


RA 


Dec 


Grating 


Frames 


Exposure 


Observed 


Quality 


PAndAS-7 


00 10 51.3 


+39 36 00.0 


g' 


6 


3600s 


2008-08-11, 2008-09-05 


0.40 - 0.60" 








i> 


5 


2400s 


2008-08-11, 2008-08-31 


0.55 - 0.65" 








R831 


6 


5400s 


2010-07-20 


< 0.75" 


PAndAS-8 


00 12 52.4 


+38 17 48.0 


g' 


6 


3600s 


2008-09-03 


0.50 - 0.55" 








%> 


5 


2400s 


2008-09-03 


0.40 - 0.50" 








R831 


6 


5400s 


2010-07-20, 2010-08-15 


< 0.75" 


Gl 


00 32 46.5 


+39 34 40.7 


R831 


9 


1350s 


2010-07-20 


< 0.75" 


MGC1 


00 50 42.5 


+32 54 58.8 


R831 


9 


1800s 


2011-08-02 


< 0.75" 



filter and five with the i' filter, arranged in a 3 x 2 dither 
pattern with a step size of 5" designed to eliminate the gaps 
between the CCDs and provide continuous coverage of our 
field. Exposure durations were 600s per image for the g' filter 
and 480s per image for the i' filter. 

We reduced our data using the CMOS software pack- 



age in iraf. Appropriate bias and flat-field images obtained 
as part of the standard CMOS baseline calibrations were 
downloaded from the Gemini science archive and then ap- 
plied to each exposure using the GIREDUCE task. The three 
CCD frames in a given exposure were next mosaicked into 
a single frame using the GMOSAIC task, and the five or six 



© 2012 RAS, MNRAS 000, [THM] 




Figure 2. Central 1' X 1' thumbnails from our combined 2400s i'-band CMOS images of PAndAS-7 (left) and PAndAS-8 (right). The 
image quality is 0.57" and 0.45", respectively. North is up and east is to the left. Longslit positioning and orientation is indicated. 



mosaicked frames for a given filter were then combined into 
a single image using the IMCOADD task. 

Fig- [U shows the central 1' x 1' regions of our combined 
2400s i'-band images of the two clusters. Stellar objects in 
the image of PA-7 have a FWHM « 7.8 pixels, or 0.57", 
while those in the image of PA-8 have a FWHM ~ 6.2 pixels, 
or 0.45". Although crowding precludes precise photometry 
in the innermost regions of both clusters, each has a suffi- 
cient number of individual members surrounding its centre 
to produce a useable colour- magnitude diagram (CMD). 

We performed photometric measurements on our com- 
bined images following a procedure identical to that which 
we previously applied to GMOS imagi ng of the exceptionall y 
remote M31 globular cluster MGC1 (|Mackev et al.ll2010bl ). 
In brief, we fit a point-spread function (PSF) model to each 
combined image us ing the stand- alone versions of daophot 
ii and ALLSTAR II |Stetsonlll987l ). For any given image we 
began by conducting a first pass of object detection with 
daophot, and selected ~ 75 relatively bright, isolated stars 
from the list of detections to construct an initial PSF. Pre- 
vious experience with MGC1, along with some additional 
experimentation on our present set of observations, showed 
that a 'penny2' PSF model (Gaussian core plus Lorentzian 
wings) varying quadratically across the GMOS field of view 
provided excellent results. We iteratively removed from the 
PSF list any stars with an error in their fit of more than 
three times the average, redefining the model each time. 

After reaching convergence (at which point there were 
still typically > 60 stars defining the PSF) we used ALLSTAR 
to subtract from the image any stars neighbouring those 
in the PSF list. We then used the now completely isolated 
PSF stars on this subtracted image to recalculate and fur- 
ther refine the PSF model. Next, we applied this model to 
the original image using ALLSTAR and subtracted all known 
stars. This subtracted image was then run back through 
DAOPHOT in order to find faint objects missed in the first 
detection pass. Finally, we took the original image and our 
final PSF model and used ALLSTAR to perform photometric 
measurements on the complete list of detected objects. 



In order to eliminate non-stellar objects and stars with 
poor photometry, we passed the resulting list of measure- 
ments through several quality niters based on parameters 
calculated by ALLSTAR. Specifically, we filtered objects by 
the x 2 °f their PSF fit, their estimated photometric er- 
ror, and their measured sharpness relative to the PSF 
model. Both clusters lie at moderately low Galactic latitude 
(b ~ —20°), resulting in plenty of foreground stars scattered 
across both GMOS fields that allowed us to define, in each 
filter, suitable ranges in these parameters as a function of 
magnitude for well- measured point sources. 

To calibrate our photometry, we decided to depart from 
the procedure we applied to our previous observations of 
MGC1. Instead of using the baseline Gemini calibration 
imaging of one standard field (containing ~ 5 stars) per 
cluster per night, we cross-matched our list of stellar de- 
tections against the PAndAS catalogue in order to derive 
suitable transformations from instrumental Gemini/GMOS 
magnitudes onto the CFHT/MegaCam photometric system. 
We excised circular regions of radius 15" about each clus- 
ter to avoid problems associated with crowding; across the 
PA-7 field we matched 131 stars while in the PA-8 field we 
matched 109 stars. This large number of stars, covering a 
wide range of colours 0.0 < (g — i)cfht < 4.0, allowed for 
a more robust transformation than we were able to achieve 
for MGC1 by employing four standard stars to calibrate our 
photometry onto the SDSS system. 

An example of our transformation procedure is shown 
in Figure [3] As noted previously, the GMOS g' and i! niters 
are similar, but not identical, to the g' and i niters em- 
ployed by CFHT. To a first approximation, we found both 
the g-band and i-band transformations to be linear. How- 
ever, some experimentation revealed a mild deviation at the 
faint end for both filters, which we accounted for by adding 
a small exponential term at faint magnitudes. The necessity 
of this additional term is not evident to the eye in Figure |3j 
but was readily identified after, for example, trying to match 
globular cluster fiducial sequences to the colour-magnitude 
diagrams for PA-7 and PA-8 (Section l3~2j) . 
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Figure 3. An example of the photometric transformation from instrumental Gemini/GMOS magnitudes to the calibrated (PAndAS) 
CFHT/MegaCam photometric system, for PA-7 i-band observations. The leftmost panel shows the relationship between i'Gemini and 
*CFHT- The transformation is nearly linear, except for a small exponential term that is required for faint detections. The solid line 
represents the linear + exponential relationship; the dotted line is the linear relationship only. The middle panels show the resulting 
residuals as a function of the instrumental i'Gemini magnitude and (g 1 — i')Gemini colour. In the upper of these panels the solid points are 
for the blue half of the detections and the open points for the red half; in the lower panel for clarity we plot only detections possessing 
uncertainties on the residuals < 0.05 mag. Both panels show that a small colour term is required - we employ a simple prescription 
as marked (dashed line). The rightmost panels show the residuals after the full transformation has been applied; these are centred on 
zero with rms of a few xO.Ol mag. In the central and rightmost panels we shade the region spanned by the cluster CMDs - although to 
avoid crowding uncertainties we excised most cluster members from the matched detection lists before deriving the transformations. 



Table 2. Coefficients for transforming our instrumental Gemini/GMOS photometry into the CFHT/MegaCam system. 



Cluster and 


Linear Term 


Exponential Term 




Colour Term 




Passband 


a 


P 


7 


5 


c 


«i 








PAndAS-7 g-band 


1.0081 


9.5523 


0.0214 


13.0 


1.1072 


0.000 


0.000 






PAndAS-7 i-band 


1.0073 


9.8127 


0.0238 


12.0 


0.2743 


0.025 


-0.050 


1.80 


3.60 


PAndAS-8 g-band 


1.0080 


9.5789 


0.0260 


13.3 


0.3417 


-0.030 


0.020 


1.00 


1.80 


PAndAS-8 i-band 


1.0001 


9.7305 


0.0338 


12.0 


0.4752 


0.025 


-0.020 


1.60 


2.80 



As in Figure O plots of the resulting residuals revealed 
that further small corrections involving (g' — £') Gemini were 
required. The implied form of these was consistent across 
the g- and i-band photometry for both clusters - a constant 
offset at blue colours, moving linearly to another constant 
offset at red colours. 

In summary, our final transformations took the form: 



mtrans = OBl' + ft + f](m') - %{ C ') 



(1) 



Here, m! is the instrumental GMOS g' or i' magnitude 
and c' is the instrumental (g' — i') colour, mt ran s is the 
final transformed magnitude on the calibrated (PAndAS) 
CFHT/MegaCam scale, a and /3 are coefficients defining 
the linear term, and ri(m') and x( c ') are t ne exponential 
and colour terms, respectively: 



ii' if m' < S 

ri(m = < 
V ; \y{m' + 8y 



X(c') 



(«2- 



K2 



Kl) 



K 3 ) + Kl 



if C < K3 

if K3 ^ c' ^ K4 

if c' > K4 . 



The coefficient 5 represents the cut-on magnitude for the ex- 
ponential term, while K3 and K4 are the blue and red colours 
of the knees in the colour term. We list the derived coeffi- 
cients for each of our photometric transformations in Table 
[2] As in Figure [3] application of these transformations re- 
sulted in rms residuals at a level of a few xO.Ol mag across 
the full range of colours in both filters. 

The MegaCam i-band filter broke in June 2007 and was 
replaced in October 2007; these two filters are not identical 
at a level of up to ±0T m ag across the full colour range (e.g. 
iMcConnachie et aLllioirl . The PAndAS imaging covering 
PA-7 and PA-8 was taken in late 2006 and our transforma- 
tions are hence to the original MegaCam i-band system. 
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2.2 Spectroscopy 

We obtained longslit spectra of PA-7 and PA-8 with GMOS- 
N on two nights in July and August 2010 with the aim of 
measuring radial velocities for these two faint clusters. The 
data were obtained in queue mode as part of program GN- 
2010B-Q-19 (PI: Mackey) during clear conditions and with 
seeing better than 0.75", as outlined in Table [T] We em- 
ployed the R831 grating operating at central wavelengths 
near 8500 A, in combination with a 0.75" slit width and 
the RG610 blocking filter, to obtain a (measured) resolution 
R » 3150 around the Can triplet at 8498 A, 8542 A, and 
8662 A. We binned the CCD by a factor two in both direc- 
tions to achieve a spatial resolution of 0.1454" per pixel and 
a dispersion of 0.68 A per pixel. PA-7 is a loosely aggregated 
cluster with no well defined core, so we placed the slit on the 
brightest central red giant member and oriented it at a po- 
sition angle of 167° east of north in order to cover the next 
brightest giant plus at least one fainter member. For PA-8 
we oriented the slit at a position angle of 28° east of north in 
order to cover the unresolved cluster core and the brightest 
resolved red giant member. Positioning of the slit for both 
targets is indicated in Figure [2] 

For each cluster we obtained six 900s exposures, split 
into two groups of three observed at central wavelengths 
of 8450 A and 8500 A respectively. This strategy served to 
fill the wavelength coverage gaps between the three GMOS 
CCDs as well as mitigating the effects of hot pixels and 
bad columns. Within each group the three exposures were 
dithered by ±15" along the slit in order to fill the gaps 
created by the stabilising bridges on the GMOS longslit 
mask, further mitigate the impact of any detector defects, 
and to help reduce systematics during the subtraction of at- 
mospheric emission lines. We bracketed each group of three 
science frames with arc-lamp exposures to ensure an accu- 
rate wavelength calibration. Flat-field images were taken for 
each central wavelength setting. 

In addition to PA-7 and PA-8 we observed the very 
bright and concentrated M31 globular clusters Gl and 
MGC1 to serve as radial velocity templates. Gl is listed 
with a velocity V r = —332 ± 3 kms" 1 in the R evised 
Bologna Catalogue (RB C, V4.0; iGalleti et alj|2004 ). while 
lAlves-Brito et all l|2009h measured K = -355 ±2 kms" 1 for 
MGC1 from a high resolution spectrum. As outlined in Table 
[T] our observations of Gl were obtained in July 2010 as part 
of the same program (GN-2010B-Q-19) under which PA-7 
and PA-8 were targeted, while our spectra of MGC1 were 
obtained in August 2011 as part of program GN-2011B-Q- 
61 (PI: Mackey). For both clusters we employed an identical 
set-up to that described above for PA-7 and PA-8, except 
that the slit was centred on the unresolved core of each tar- 
get. The brightness of the two clusters necessitated shorter 
individual exposure durations of 150s for Gl and 200s for 
MGC1. We obtained nine frames for each target, split into 
three groups of three observed at central wavelengths of 
8450 A, 8500 A, and 8550 A. As before, within each of these 
groups the three exposures were dithered by ±15" along the 
slit, and each group was bracketed with arc-lamp exposures. 

As with the imaging data, we reduced our spectro- 
scopic observations using the GMOS package in iraf. For a 
given cluster, we first created normalised flat-fields for each 
of the central wavelength settings using the GSFLAT task. 




8450 8500 8550 8600 8650 8700 
Wavelength (A) 



Figure 4. Combined normalized spectra around the Call triplet 
for PA-7 and PA-8, bracketed by our two velocity template clus- 
ters Gl and MGC1. The vertical dotted lines indicate the posi- 
tions of the Can lines in the spectrum for Gl. Note that PA-7 
and PA-8 are clearly blue-shifted with respect to these lines. 



These normalised flat-fields were then applied to the science 
frames, along with an appropriate bias image from the Gem- 
ini archive, using the GSREDUCE task. We also ran gsreduce 
on the arc-lamp exposures, but did not apply flat-fields to 
these frames. Next, we ran GSWAVELENGTH on the reduced 
arcs in order to establish a wavelength calibration at all 
points along the slit for any given science frame, and then 
applied these calibrations and rectified the science frames 
with GSTRANSFORM. We carefully examined the positions of 
the sky lines on each frame to check that the calibration 
was accurate, particularly near the object spectrum and the 
Call triplet. We applied GSSKYSUB to perform a sky subtrac- 
tion on each of the calibrated science exposures, and then 
extracted ID object spectra from the sky-subtracted frames 
using the GSEXTRACT task. This step was straightforward for 
the two reference clusters - for each one we simply extracted 
the spectrum of the bright unresolved core. For PA-8 we ex- 
tracted the spectrum of the unresolved core as well as that 
of the bright red giant member we also placed on the slit. 
For PA-7 we attempted to extract spectra for all three of the 
red giants we placed on the slit; however only the brightest 
had sufficient signal to provide a useable spectrum. 

Finally, we corrected each extracted spectrum to the he- 
liocentric frame using the iraf tasks RVCORRECT and dop- 
COR, and then median-combined the individual exposures of 
a given target into a final spectrum using the SCOMBINE task. 
Best results, including cosmic-ray removal, were achieved by 
scaling the input spectra to a common flux level and then, 
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when combining them, weighting each by its median flux 
value and applying a sigma-clipping rejection algorithm. 

Figure [3] shows our combined spectra in the region of 
the Call triplet for PA-7, PA-8 and the two reference clus- 
ters; each has been normalized using the iraf task con- 
tinuum. We estimate that in this part of the spectrum the 
signal-to-noise is w 7 per resolution element for PA-7, and 
« 20 per resolution element for PA-8. Even so, the Can 
lines are clearly visible in both cases - their positions sug- 
gest that both PA-7 and PA-8 have larger blue-shifts than 
either template cluster. 



velocity dispersion to decrease with projected radius R as: 



a v (R) = 152-0.90 



R 

kpc 



kms 



(2) 



3 RESULTS AND ANALYSIS 



3.1 Radial velocities 



We used the iraf task FXCOR to determine radial velocity 
estimates for PA-7 and PA-8 by cross-correlating their spec- 
tra again st those of the temp late clusters Gl and MGC1 
(see e.g.. iTonrv fc Davidll979h . Because of the presence of 
strong telluric absorption features at several different places 
within our spectral coverage, as well as the presence of non- 
zero residuals from the sky subtraction, we chose to cross- 
correlate just the region covering the Call triplet with a 
~ 50 A buffer at either end (i.e., ~ 8450 - 8700 A). This 
worked well for PA-8, producing a strong isolated cross- 
correlation peak; however for PA-7 the results were poor. 
We traced the problem for PA-7 to the Ca II line at 8662 A 
which falls in a crowded region of sky emission. Due to the 
low signal-to-noise of the source spectrum, residuals from 
the sky subtraction skewed the shape (and position) of this 
Call line. By restricting the cross-correlation region to the 
Call lines at 8498 A and 8542 A (i.e., adopting the range 
~ 8450 — 8570 A) we circumvented this problem. We ver- 
ified that using just this region for PA-8, we were able to 
reproduce the velocity obtained from using all three Call 
lines to better than the uncertainty on that measurement. 

Our results are as follows: for PA-7 we obtained a helio- 
centric radial velocity V r = —433 ± 8 kms -1 , where the un- 
certainty includes the quoted errors on the velocities of the 
reference targets. The height of the cross-c orrelation peak, 
agains t either template, was « 0.5, while the lTonrv fc Davisl 
l|l979h R-parameter « 30 - both indicating a sound match 
between the target and reference spectra. For PA-8 we found 
V r = — 411±4 km s _1 , with a cross-correlation peak of height 
« 0.8, and R ~ 49 - indicating a high quality match between 
the spectra. 

These two velocities fall quite close together and are 
rather different from the systemic M31 velocity of —301 
kms -1 , suggesting that PA-7 and PA-8 are likely to be as- 
sociated with the underlying South- West Cloud rather than 
being two unrelated members of the M31 halo. Ultimate 
verification of this association will require velocity measure- 
ments for stellar members of the tidal substructure; for now 
we use a simple statistical test to quantify the significance 
of our result. 

By measuring velocities for a large number of stars 
at a variety of diff erent positions in the M31 halo, 
IChapman et al.l (|2006h estimated the (1-dimensional) halo 



Extrapolating a small amount beyond their outermost field 
at R w 70 kpc suggests that the velocity dispersion of the 
M31 halo at radii comparable to PA-7 and PA-8 (R ~ 87 
kpc) should be a v (R) ~ 70 kms -1 . Under the assumption 
that remote M31 globular clusters share the global kinemat- 
ics of the metal-poor field halo, we ask how likely it is that 
two independent objects, not members of any accreted sub- 
structure, should have velocities at least as far away from 
the M31 systemic velocity as observed for both PA-7 and 
PA-8, but within a similar margin of each other. We used 
a simple Monte Carlo model to address this question - we 
generated a very large number (> 10 6 ) of globular cluster 
pairs with velocities randomly drawn from a Gaussian distri- 
bution of width 70 kms -1 centred on Vm3i = —301 kms -1 , 
and counted the number of times that (i) both clusters had 
velocities \V ri 2 — Vm3i| *5 100 kms -1 , and (ii) the differ- 
ence between them \V ri — V r2 \ ^ 30 kms -1 . We found this 
to be a very unusual configuration, occuring just ~ 1.5% of 
the time. Note that the bounds adopted here are somewhat 
conservative to allow for the uncertainties in our radial ve- 
locity measurements, meaning that this probability is best 
treated as an upper limit. 

This simple test indicates that PA-7 and PA-8 are al- 
most certainly not unrelated members of the M31 halo, con- 
sistent with our observation that these two clusters are pro- 
jected directly on top of the South- West Cloud and strongly 
suggestive that they are both associated with this substruc- 
ture. 



3.2 Metallicity and age constraints 

Colour-magnitude diagrams for PA-7 and PA-8 are displayed 
in Figure [5] Our photometry reaches to g ~ 27 and i « 26, 
just over one magnitude below the typical horizontal branch 
level for globular clusters at the M31 distance and redden- 
ing. The CMDs are sufficiently well populated to discern 
the locus of the red giant branch (RGB) and the position 
and morphology of the HB. For both clusters the latter fea- 
ture is particularly striking, being very short and exclusively 
red - in fact, more of a red clump than a typical globu- 
lar cluster HB. In this respect the CMDs bear considerable 
resemblence to those for many of the young halo globu- 
lar clusters found in the outskirts of the Milky Way, such 
as Palomar 3, 4, and 14, Pyxis , Eridanus, and AM-1 (e.g., 
Stetson et al. Ill999l: lHilkerll2006l: iDotter. Saraiedini fc Yand 



2008al : lDotter. Saraiedini fc Andersonll201ll y These part 
lar objects possess ages several Gyr you nger than the oldest 
members of the Galact ic system (e.g., iMarin-Franch et al.l 
I2009l ; [l5otter et al"1l201Cf ). Given this, it is quite possible that 
both PA-7 and PA-8 are also young globular clusters, and 
we assess the likelihood of this below. 

The first point to note is that crowding and incomplete- 
ness do not affect the observed HB morphology for PA-7 and 
PA-8. That is, we do not believe that there are blue HB stars 
in these clusters that were preferentially missed during the 
photometric measurements. A by-eye inspection of the de- 
tected objects indicates that all stars at this brightness level 
in the vicinity of the clusters were photometred. A few were 
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Figure 5. Colour-magnitude diagrams for PA-7 and PA-8. Typical photometric uncertainties are as indicated. The upper RGB stars 
circled in red are those for which we obtained useful longslit spectra (note that for PA-8 we also observed the unresolved cluster core). 



discarded due to crowding issues but none of these were sig- 
nificantly bluer than the HB stars with good measurements. 
As further evidence of this point, we direct the interested 
reader to our observa tions of the M31 halo cluster MGC1 
|Mackev et all l2010tj ). which we took with CMOS under 
very similar atmospheric conditions. Crowding and incom- 
pleteness is more of a problem in this luminous compact 
cluster, yet we still clearly detect stars all along its blue HB. 

We estimated the metallicity, foreground reddening, 
and distance of PA-7 and PA-8 by overplotting fiducial 
sequences from several Galactic globular clusters on their 
CMDs. For a given comparison cluster we measured the ver- 
tical shift required to align the HB level with that of either 
PA-7 or PA-8, along with the horizontal shift required to 
align the colour of the red giant branches (RGBs) at the 
HB level. This horizontal shift gives the difference in colour 
excess between the two clusters; when combined with the 
vertical shift this then yields the difference in distance mod- 
ulus. Since the shape of the upper RGB is strongly sensitive 
to [Fe/H], the fiducial sequence most closely tracing this fea- 
ture indicates the cluster metallicity. 

Because the colour of the RGB at HB level is sensi- 
tive to age as well as [Fe/H] (and also, indeed, to [a/Fe]) 
we attempted, as far as possible, to adopt fiducial se- 
que nces for Galactic clusters posse s sing r ed HB morphol- 
ogy. Idem. Vanden Berg fc Stetson! l|2008h provide fiducial 
sequences for 4 globular clusters and one open clu s ter in 
the CFHT/MegaCam filter system, while lAn etafl [j2008h 
provide sequences for 17 globular clusters and 3 open clus- 
ters in the SDSS filter system. Since our final photometry 
for PA-7 and PA-8 is calibrate d to the Me gaCam system, 
ideally we would have used the Idem et al.l |2008l ) clusters 
as our templates. Unfortunately, however, a small amount 
of testing revealed that of their sample, only M71 was in 
any way close to PA-7 and PA-8 in terms of [Fe/H] and HB 
morphology. 

The sample of lAn et al.l {2008), on the other hand, pos- 
sesses a number of red HB clusters with [Fe/H] comparable 



to that for PA-7 and PA-8. Of particular relevance here are 
the outer halo objects Palomar 3 and 4. Since the SDSS pho- 
tometric system is slightly different to the CFHT/MegaCam 
sys tem, it wa s neces sary to apply suitable transformations to 
the I An et al.l |2008l ) fiducials before comparing them to our 
CMDs. We used the appropriate equations from the online 
documentation describing the MegaPipe pipeline reduction 
process for MegaCam imaging: 

<?mc = #sdss — 0.153 ((?sdss — )"sdss) 

«MC = «SDSS — 0.085 (rsDSS — «SDSS) (3) 

We also desir ed a more metal -rich template cluster than Pal 
3 or 4. In the I An et all l|200Sft sample, only M71 fits the bill; 
however employing this cluster as a template was a com- 
plex procedure. M71 is relatively nearby, meaning that its 
more lu minous RGB me mbers are saturated in SDSS im- 
ages. The[A n et alj j2008) fiducial is consequently truncated 
at a level just brighter than the HB. To circumvent this is - 
sue we adopted the M71 sequence from lClem et al.l £2008) . 
Unfortunately, these authors do not provide HB fiducials 
for their clusters, or a public archive of their photometry. 
Therefore, we had to merge the An et al. HB photometry 
onto the Clem et al. RGB fiducial. In princi ple this simply 
require s the straightforward use of Eq.[3]as in lMackev et al.l 
|2010bf ); however An et al., when comparing their photom- 
etry to that of Clem et al., found a significant disagreement 
for M71 (note that all other clusters matched to a level bet- 
ter than ~ 2%). They ascribe this discrepancy to substantial 
uncertainties in their zero points for the M71 SDSS imaging. 
To correct for this, we tested a grid of zero-point offsets in 
gri. Each combination of offsets was applied to the full An 
et al. M71 fiducial, which was then transformed using Eq. 
[31 We located the set of offsets producing the best match 
between this transformed sequence and that of Clem et al., 
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Figure 6. Galactic globular cluster fiducial sequences overplotted on the CMDs for PAndAS-7 and PAndAS-8. These have been aligned 
using the level of the HB and the colour of the RGB at the HB level, as described in the text. The upper panels show results for three 
fiducial sequences of different metallicity: Palomar 3 with [Fe/H] = —1.6, Palomar 4 with [Fe/H] = —1.35, and M71 with [Fe/H] = —0.8. 
The lower panels show the fiducial sequence and HB for Palomar 4 only, which provides a good fit to the PA-7 and PA-8 CMDs. Note 
that in both cases the Pal 4 HB extends further to the blue and has a bluer median colour than the M31 cluster. 



and then applied this to the An et al. HB photometry to 
properly transform it onto the Clem et al. M71 scale. 

The results of our template alignment are presented in 
the upper panels of Figure [6] Based on the shape of the 
upper RGB, Palomar 3 is clearly more metal-poor than 
our M31 objects, while M71 is cons i derab ly more metal- 
rich. The 2011 update of the iHarrij \ 19961 ) globular clus- 
ter catalogue lists [Fe/H] ~ —1.6 and —0.8 for these two 
objects, respectively. Palomar 4, on the other hand, pro- 
vides an excellent upper RGB fit for both PA-7 and PA- 
8, indicating that they both share a similar metallicity to 
this globular cluster. The Harris catalogue lists [Fe/H] « 

— 1.4 for Pal 4, commensurate with the value of [Fe/H] = 

— 1.41 ±0.04 (statistical) ±0.17 (systematic) measure d from 
high resolution spectroscopy bv lKoch fc Cotel |2010l ). How- 
ever, other publications suggest it may be slightly more 



m etal-rich than this - for example , the c alcium triplet work 
of lArmandroff . Da Costa fc Zinnl l|l992f) which returned a 
value [Fe/H] = —1.28 ± 0.20, or photometri c studies from 
HST imaging which find [Fe/H] « -1.3 (e.g.. IStetson et all 
ll999l ; lDotter et al.ll201ll ). Here we adopt [Fe/H] = -1.35 for 
Pal 4, and hence PA-7 and PA-8, but recognize a likely sys- 
tematic uncertainty of ±0.1 dex in this quantity. Based on 
the quality of our photometry and the separation between 
the fiducial sequences of our template clusters, we estimate 
a statistical uncertainty of ±0.15 dex on our best-fit metal- 
licity for PA-7 and PA-8. 

The Harris catalogue lists a small colour excess E(B — 
V) = 0.01 for Pal 4, and an intrinsic distance modulus (m — 
M)o = 20.18. Based on the vertical and horizontal shifts 
required to align the template fiducial with our photometry, 
the implied colour excess and distance modulus for PA-7 are 
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Figure 7. Age as a function o f metallicity and two different HB morphology metrics for a sample of 70 Galactic globular clusters, 
after iDotter et al. | <20ld . l201ll) . The location of PA-7 and PA-8 is indicated with a black point. 



E(B - V) = 0.11 ± 0.02 and (m - M) = 24.47 ± 0.07. For 
PA-8 we find E(B - V) = 0.12 ± 0.02 and (m - M) = 
24.58 ± 0.07. Note that in deriving these q uantities we have 
emplo yed the extinction coefficients from IStoughton et al.l 
(2002), specifically that Ai = 2.086 E(B - V). Our quoted 
uncertainties reflect only the precision with which we are 
able to align the fiducial sequences with our photometry - 
i.e., the uncertainties are obtained directly from those asso- 
ciated with the measured horizontal and vertical shifts. 

Our measured colour excesses are in good agree- 
ment with predictions from th e reddening maps of 
ISchlegel, Finkbeiner fc Davis! i|l99Sft : E(B -V) = 0.09 and 
0.11 for PA-7 and PA-8, respectively. Similarly, our dis- 
tance moduli match quit e well the canonica l M31 value 
(m - M) = 24.47 (e.g., iMcConnachie et al.ll2005T ). It is 
possible, at marginal significance, that PA-8 could be up to 
~ 0.1 mag more distant than PA-7. This would imply a sig- 
nificant depth of up to ~ 40 kpc to the South- West Cloud 
along the line of sight. However, we must bear in mind the 
lcr random uncertainties of 0.07 mag on each of the distance 
measurements, and that these do not include any systematic 
components due to, for example, the photometric calibra- 
tions, transformation of the Pal 4 fiducial to the MegaCam 
photometric system, etc. It is therefore safe to say the dis- 
tance estimates are consistent with each other, although it 
would not be surprising if there were a mild line-of-sight 
extension to the South- West Cloud. 

The distinctive HB morphology observed for our two 
M31 clusters provides an upper limit to t heir ages. To 
demonstrat e this we consider the recent work of lDotter et al.l 
(|20ldl201ll ) who, using HST/ACS photometry, investigated 
HB morphology in nearly half the known Galactic globular 
clusters as a function of their other properties (in particu- 



lar their ages and metal abundances). We rep roduce one of 
the k ey results of this work (Figure 11 from IDotter et all 
2011) in Figure [7] This shows the relationship between a 
cluster's age, HB morphology and [Fe/H]. HB morphology 
is paramete rized using two diff erent metrics. The first, in- 
troduced bv lDotter et al.1 i|2010T ), is the difference A(V - I) 
between the median colour of stars on the HB and the 
median colour of stars on the RGB at the HB level. The 
secon d is the commonly used rati o (B — R)/(B + V + R) 
(e.g.. lLee. Demarque fc Zinnlll994h . where B (alternatively, 
R) represents the number of HB stars lying to the blue (red) 
of the instability strip, and V the number of RR Lyrae vari- 
ables (i.e., the number of HB stars lying on the instability 
strip). 

Irrespective of which HB metric is considered, it is clear 
that in the metallicity range —1.6 < [Fe/H] < —1.0 the mor- 
phology of the HB correlates strongly with the cluster age. 
The oldest objects in the Galactic system have very blue 
horizontal branches; those several Gyr younger have very 
red horizontal branches; and there is a relatively smooth 
transition in between. This metallicity range is the "sweet 
spot" for sensitivity to age - in more metal-rich clusters 
the HB morphology is driven predominantly by [Fe/H] such 
that it is exclusively red irrespective of age, while in more 
metal-poor clusters the correspondence with age is less clear 
as only very old systems are observed with [Fe/H] < —1.6. 
Note that while there is some evidence that additional pa- 
rameters - such as the central density of a cluster, or its 
helium abundance - may also affect HB morphology, in gen- 
eral these appear to be c onsiderably less in fluential than 
both [Fe/H] and age (e.g.. IDotter et al.ll201ol ). 

At [Fe/H] « -1.35, our two M31 clusters fall directly 
in the region where the HB is most strongly sensitive to 
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cluster age. Their very red HB morphologies indicate that 
they are likely to be at least ~ 2 — 3 Gyr younger than the 
oldest Galactic globulars. To demonstrate this explicitly we 
derive the parameters A(V - I) and (B - R)/(B ± V ± R) 
for both objects. The latter metric is straightforward - all 
HB stars are red so the HB-type converges to —1.0. Al- 
though we do not explicitly know the position of the insta- 
bility stnrj_oiiJhe_^fr2 : JheHB ; type of Pal 4 is also -1.0 
(e.g., iMackev fc van den Berghll2005l ) and the lower panels 
of Figure [6] show that all HB stars in our two M31 clusters 
are redder than the bluest stars on the Pal 4 HB. 

To calculate A(V — 7) for a given cluster, we drew a 
box around the HB and calculated the median colour of 
the stars within it. Note that because we observed only a 
small number of HB stars (« 10) in each system, the me- 
dian colour is somewhat sensitive, at a level of roughly ±0.03 
mag, to the precise boundaries of the isolating box. In prin- 
ciple we should have followed a similar procedure to deter- 
mine the median colour of the RGB at the HB level; however 
we possessed too few RGB stars for this calculation to be 
reliable. Instead, we adopted the colour of the best-fitting 
(shifted) Pal 4 fiducial at the HB level, and conservatively as- 
sumed a total uncertainty of ±0.05 mag on the difference be- 
tween this quantity and the median colour of the HB. These 
two measurements gave us the quantity A(g — i), which we 
converted t o A(V — I) u s ing th e transfor mation procedure 
outlined bylMartin et all (120061 ) (see also llbata et al.ll2007l; 
iHuxor et alj|2008j ). For~PA-7 we measured A(V - I) = 0.14 
and for PA-8 we found A(V - I) = 0.12. Both values 
are smaller than t hat seen for Pal 4, A(V — I) = 0.183 
i|Dotter et al.ll2010l ). consistent with the lower panels in Fig- 
ure [6] 

Plotting the two HB metrics on Figure [7] confirms that 
PA-7 and PA-8 lie in a region inhabited only by Galactic 
globular clusters with ages at least 2 Gyr younger than those 
for the oldest members of the Milky Way system. 

While this argument provides an upper bound to the 
ages of PA-7 and PA-8, it says nothing about how young 
these objects could be. To investigate this we fi t Dartmouth 
isochro nes and synthetic horizontal branches (|Dotter et al.l 
l2008bl ) of different ages to our CMDs. We assume [Fe/H] = 
— 1.35 and an a-element abundance equivale nt to that ob- 
served for Palomar 4 - i.e., [a/Fe] » ±0.4 (|Koch fc Cotel 
l2010h . Note that our results are rather insensitive to the 
precise value adopted for the latter parameter. 

Based on the level of the isochrone sub-giant branches, 
our faint detection limits of (g, i) ~ (27.0, 26.0) formally 
correspond to a minimum age of w 2 Gyr for PA-7 and 
PA-8. However, the implied distance moduli and foreground 
reddening levels necessary to correctly align such young 
isochrones on the cluster CMDs suggest that their ages are, 
in all probability, significantly greater than this. Figure [8] 
explicitly demonstrates this result. We plot isochrones rang- 
ing in age from 2.0 Gyr to 12.5 Gyr, shifted according to the 
best-fit distance moduli and foreground extinction values de- 
rived above from aligning the Palomar 4 fiducial. Because 
the colour of the RGB at the HB level is mildly sensitive 
to age in the sense that younger clusters have somewhat 
bluer RGB positions, applying these shifts does not lead 
to the youngest isochrones (or indeed their corresponding 
HB models) being well aligned; significantly larger values of 
E(B—V) and (m— M)q would be necessary to correct this. In 



numerical terms, aligning the 2.0 Gyr isochrone on the PA-7 
CMD implies E(B - V) = 0.16 and (m - M) = 24.7, and 
for PA-8 E(B-V) = 0.17 and (m-M)o = 24.8. These fore- 
ground reddenin g levels are significan tly greater than those 
predicted by the lSchlegel et aL (119981 ) maps; note, however, 
that such a discrepan cy is not out of the question (see e.g., 
IMackev et al.l l2010bl ) . More problematic would be the co- 
incidence of this discrepancy with the unusually large im- 
plied distance moduli, which would place the clusters, re- 
spectively, « 90 and 130 kpc behind the central disk of 
M31. Arguably a more satisfactory resolution is achieved 
by adopting an age in the range « 6 — 10 Gyr, in which case 
both E(B — V) and (m — M)o are much more consistent 
with expected values. 

Isochrones with age 12.5 Gyr can be aligned nicely with 
the CMDs if, in both cases, E(B - V) is about 0.02 mag 
smaller and (m — M)q about 0.1 mag smaller than our previ- 
ous best-fit values. These are certainly within the bounds of 
acceptability; however as discussed previously the observed 
HB morphologies of our two M31 clusters exclude such an 
old age (Figure [7} . 



4 DISCUSSION 

The measurements presented here for PA-7 and PA-8 rep- 
resent a new step in improving our understanding of the 
assembly of the outer halos and globular cluster systems of 
M31 and the Milky Way. As discussed in Section T3. II the ve- 
locities of these objects are sufficiently similar to each other, 
and, given their large projected galactocentric radii, suffi- 
ciently different from the M31 systemic velocity, that it is 
very likely that they are not independent clusters but rather 
physically associated with the underlying South- West Cloud 
(although with the caveat that to establish the link directly 
will r equire velocit i es for stellar members of this overden- 
sity). ICollins et al.l (2009) have previously shown that the 
extended M31 cluster HEC12 (also sometimes called EC4) 
shares a common velocity with the underlying halo feature 
known as Stream C. Apart from this, our measurements for 
PA-7 and PA-8 represent the strongest confirmation to date 
that some remote M31 halo globular clusters are members of 
stellar substructures - despite this commonly being assumed 
to be the case as a result of the statistically significant cor- 
relation between the spatial positions of man y such objects 
and a number of the field halo substructures l|Mackev et al.l 
l2010al ). 

In this context the resolved properties of PA-7 and PA- 
8 take on added significance. We have shown that both have 
[Fe/H] = — 1.35±0.15, and both possess a strikingly red hor- 
izontal branch morphology for this metallicity that implies 
their ages must be at least ~ 2 Gyr younger than the oldest 
Galactic globular clusters. With such characteristics PA-7 
and PA-8 would be unambiguously classified a s members of 
the "y oung halo" subsystem in the Milky Way. iPerina et al.l 
( 2012) have recently demonstrated that many halo globular 
clusters in M31 exhibit an even stronger second parame- 
ter effect than do the young halo clusters in the Milky Way. 
They attribute this to the M31 objects having formed ~ 1 — 2 
Gyr later than their Milky Way counterparts, which is en- 
tirely consistent with our conclusions for PA-7 and PA-8. 

As discussed in the introduction to this paper, there is a 
wealth of circumstantial evidence that the young halo pop- 
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Figure 8. Dartmouth isochrones and synthetic HBs overplotted on the CMDs for PA-7 and PA-8. The models have [Fe/H] = —1.35 
and [ct/Fe] = +0.4, and a variety of ages as indicated. For illustrative purposes (see text) they have been shifted according to the 
foreground extinction levels and distance moduli derived from aligning the Palomar 4 fiducial. For clarity we do not plot a synthetic 
HB for the 10 — 12.5 Gyr isochrones; this falls somewhat fainter (< 0.1 mag) than the cluster HB levels. 



ulation of globular clusters has been accreted into the outer 
Galaxy from now-defunct satellite systems; however, direct 
observations of this process in action are limited to a handful 
of members and ex-members of the Sagittarius dwarf. Our 
measurements provide new, independent evidence for close 
analogues of Galactic young halo clusters being accreted into 
the halo of a large spiral galaxy, in a manner entirely consis- 
tent with our understanding of the origin of the outer parts 
of the Milky Way globular cluster system. 

The complementary viewpoint is also interesting. If we 
accept a priori that the available evidence renders the ac- 
cretion picture largely correct for the Milky Way (despite 
very few clusters having been directly associated with stellar 
streams), then our observations of PA-7 and PA-8 are signif- 
icant in that they explicitly demonstrate that objects with 
very similar, if not identical, properties to the accreted mem- 
bers of the Galactic system are also b eing accreted into M31 . 
In f act, combined with the results of IMackev et al.l (|2010al ) 
and IPerina et all (|2012h . our results suggest that this pro- 
cess has occurred on a grand scale in M31. This would add 
credence to the idea that such processes play a key role in 
determining the composition of globular cluster systems in 
large spiral galaxies in general. It is of considerable interest 
to extend our study to include resolved observations of a 
more statistically robust number of globular clusters lying 
on (and off) a wide variety of M31 substructures in order to 
test how far the similarities with the outer Galactic globular 
cluster system stretch. 

A number of autho r s (e.g., Marin-Franch et alj 120091 ; 
iPorbes fc Bridges! l20ld : iDotter et all 1201(1 l201ll ) have 
demonstrated that the age-metallicity relation (AMR) for 
Galactic globular clusters exhibits two distinct loci - one at 
almost constant old age (~ 13 — 13.5 Gyr) and a second 
branching to younger ages at [Fe/H] > —1.5. The younger 
branch is composed mainly of young halo clusters includ- 



ing several of those associated with the Sagittarius dwarf. 
It is interesting to consider PA-7 and PA-8 in this setting. 
Our measured [Fe/H] « —1.35 places these two objects 
just above the branch-point on the Galactic globular cluster 
AMR. At this metallicity, Galactic globular clusters are seen 
with ages up to ~ 3.5 Gyr younger than the oldest members 
of the system (see e.g., Figure [7]). As discussed in Section 
3.21 we formally constrain the ages of PA-7 and PA-8 to 
lie in the range ~ 2.0 — 11.0 Gyr; however in all likelihood 
they have ages older than « 6 Gyr. Thus it is quite possi- 
ble that these two objects fit well with the younger branch 
of the Galactic globular cluster AMR, in which case their 
galaxy of origin - the progenitor of the South- West Cloud 
- should have had a cluster enrichment history similar to 
those satellite galaxies which have deposited their members 
in the outer Milky Way. 



On the other hand, it is also possible that PA-7 and 
PA-8 have younger ages than any known Galactic globular 
cluster at their metallicity, in which case the implied clus- 
ter enrichment history of the SW Cloud progenitor would 
be comparatively rather slow, perhaps akin to that seen in 
the SMC. Note that a handful of globular clusters with ages 
as young as w 6 — 9 Gyr are seen in the Milky Way (e.g., 
Palomar 1 and 12, and Terzan 7); however these objects are 
all considerably more metal-rich than PA-7 and PA-8, with 
-0.9 < [Fe/H] < -0.5. At least two (Pal 12 and Ter 7) are 
associated with the Sagittarius dwarf. To more accurately 
constrain the ages of PA-7 and PA-8, and thus learn some- 
thing about the cluster AMR of their host galaxy, would 
require observations approaching their main-sequence turn- 
offs on the CMD - which at present can only be achieved 
with HST. Failing this, multicolour broadband integrated 
photometry could be used to place less stringent constraints 
on the ages of these two objects. 

© 2012 RAS, MNRAS 000,[THT4] 
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5 SUMMARY & CONCLUSIONS 

In this paper we have presented spectroscopy and colour- 
magnitude diagrams for two newly-discovered globular clus- 
ters in the remote halo of M31. These two objects, PA-7 and 
PA- 8, project onto a field substructure known as the South- 
West Cloud, leading previous authors (e.g., iMackev et al.l 
2010a) to suggest that they have probably been accreted 
into the M3I system. We measure radial velocities of V r = 
-433 ± 8 kms" 1 for PA-7 and -411 ± 4 kms" 1 for PA-8. 
These are rather similar to each other, especially given their 
moderately large separation from the M3I systemic velocity. 
Based on the known velocity dispersion of metal-poor stellar 
populations in the outer parts of the M3I halo, we used a 
simple Monte Carlo model to show that observing velocities 
similar to those for PA-7 and PA-8 is rather unlikely if these 
objects are completely independent, with a probability of 
< 1.5%. This result reinforces the link between PA-7, PA-8 
and the SW Cloud and provides strong evidence that these 
clusters have indeed been accreted into the M31 halo. 

Our colour-magnitude diagrams indicate that both PA- 
7 and PA-8 have metallicities [Fe/H] = -1.35 ± 0.15. A 
reference fiducial for the remote Galactic globular cluster 
Palomar 4 provides an excellent fit to the red giant branches 
in both clusters. From the horizontal and vertical shifts nec- 
essary to correctly align the Pal 4 sequence we derive fore- 
ground extinctions consistent wi t h tho se predicted by the 
reddening maps of lSchlegel et al.l 1 19981 ) and distance mod- 
uli commensurate with the canonical M31 distance. PA-8 
is up to ~ 0.1 mag (~ 40 kpc) more distant than PA-7 
at marginal significance; if this observation is correct then it 
implies a moderate line-of-sight depth to the SW Cloud that 
could, in future, be used to help constrain the geometry of 
the orbit of its progenitor galaxy. 

The most striking aspect of the PA-7 and PA-8 colour- 
magnitude diagrams is the extremely short, red horizontal 
branch morphology exhibited by both clusters. In the Galac- 
tic halo, clusters with similar [Fe/H] and HB morphology are 
exclusively > 2 Gyr younger than the oldest members of the 
system. This implies that PA-7 and PA-8 are very likely to 
have similarly young ages. Formally, our lower age limit for 
both clusters is « 2 Gyr; however we provide circumstantial 
arguments that they are probably > 6 Gyr in age. Our obser- 
vations provide strong evidence for young globular clusters 
being accreted into the halo of a large spiral galaxy, in a 
manner entirely consistent with our picture for the assem- 
bly of the outer Milky Way globular cluster system. These 
results further add credence to the idea that such processes 
play a central role in determining the composition of globu- 
lar cluster systems in large spiral galaxies in general. 

ACKNOWLEDGMENTS 

We would like to thank Aaron Dotter for his help in con- 
structing Figure (7) We appreciate suggestions from the ref- 
eree which helped improve the manuscript. ADM is grateful 
for support by an Australian Research Fellowship (Grant 
DP1093431) from the Australian Research Council. AMNF, 
ADM, and APH acknowledge support by a Marie Curie Ex- 
cellence Grant from the European Commission under con- 
tract MCEXT-CT-2005-025869 during which this work was 
initiated. GFL thanks the Australian Research Council for 



support through his Future Fellowship (FT100100268) and 
Discovery Project (DP110100678). 

This paper is based on observations obtained at the 
Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooper- 
ative agreement with the NSF on behalf of the Gemini part- 
nership: the National Science Foundation (United States), 
the Science and Technology Facilities Council (United King- 
dom), the National Research Council (Canada), CONICYT 
(Chile), the Australian Research Council (Australia), Min- 
isterio da Ciencia e Tecnologia (Brazil) and SECYT (Ar- 
gentina). These observations were obtained under programs 
GN-2008B-Q-22, GN-2010B-Q-19, and GN-2011B-Q-61. 

REFERENCES 

Alves-Brito A., Forbes D.A., Mendel J.T., Hau G.K.T., 

Murphy M.T., 2009, MNRAS, 395, L34 
An D., et aL, 2008, ApJS, 179, 326 

Armandroff T.E., Da Costa G.S., Zinn R., 1992, AJ, 104, 
164 

Bellazzini M., Ferraro F.R., Ibata R.A., 2003, AJ, 125, 188 
Chapman S.C., Ibata R., Lewis G.F., Ferguson A.M.N. , 

Irwin M., McConnachie A., Tanvir N., 2006, ApJ, 653, 

255 

Clem J.L., Vanden Berg DA., Stetson P.B., 2008, AJ, 135, 
682 

Cohen J.G., 2004, AJ, 127, 1545 

Collins M.L.M., et al., 2009, MNRAS, 396, 1619 

Da Costa G.S., Armandroff T.E., 1995, AJ, 109, 2533 

Dotter A., Sarajedini A., Yang S., 2008a, AJ, 136, 1407 

Dotter A., Chaboyer B., Jevremovic D., Rostov V., Baron 

E., Ferguson J.W., 2008b, ApJS, 178, 89 
Dotter A., et al., 2010, ApJ, 708, 698 
Dotter A., Sarajedini A., Anderson J., 2011, 738, 74 
Ferguson A.M.N., Irwin M.J., Ibata R.A., Lewis OF., Tan- 
vir N.R., 2002, AJ, 124, 1452 
Forbes DA., Bridges T., 2010, MNRAS, 404, 1203 
Galleti S., Federici L., Bellazzini M., Fusi Pecci F., Macrina 

S., 2004, A&A, 416, 917 
Harris W.E., 1996, AJ, 112, 1487 
Hilker M., 2006, A&A, 448, 171 

Hook L, J0rgensen I., Allington-Smith J.R., Davies R.L., 
Metcalfe N., Murowinski R.G., Crampton D., 2004, PASP, 
116, 425 

Huxor A.P., Tanvir N.R., Irwin M.J., Ibata R., Collett J.L., 
Ferguson A.M.N. , Bridges T., Lewis G.F. 2005, MNRAS, 
360, 1007 

Huxor A. P., Tanvir N.R., Ferguson A.M.N. , Irwin M.J., 
Ibata R.A., Bridges T., Lewis G.F., 2008, MNRAS, 385, 
1989 

Huxor A.P, et al., 2011, MNRAS, 414, 770 
Ibata R.A., Gilmore G., Irwin M.J., 1994, Nature, 370, 194 
Ibata R.A., Gilmore G., Irwin M.J., 1995, MNRAS, 277, 
781 

Ibata R., Martin N.F., Irwin M., Chapman S., Ferguson 
A.M.N., Lewis G.F., McConnachie A.W., 2007, ApJ, 671, 
1591 

Keller S.C., Mackey A.D., Da Costa G.S., 2012, ApJ, 744, 
57 

Koch A., Cote P., 2010, A&A, 517, A59 

Law D.R., Majewski S.R., 2010, ApJ, 718, 1128 



© 2012 RAS, MNRAS 000, HE! 



14 A. D. Mackey et al. 



Lee Y.-W., Demarque P., Zinn R.J., 1994, ApJ, 423, 248 
Lewis, G.F., et al., 2012, submitted 
Mackey A.D., Gilmore G.F., 2004, MNRAS, 355, 504 
Mackey A.D., van den Bergh S., 2005, MNRAS, 360, 631 
Mackey A.D., et al., 2010a, ApJ, 717, Lll 
Mackey A.D., et al., 2010b, MNRAS, 401, 533 
Marin-Franch A., et al., 2009, ApJ, 694, 1498 
Martin N.F., Ibata R.A., Irwin M.J., Chapman S., Lewis 
G.F., Ferguson A.M.N. , Tanvir N., McConnachie A.W., 
2006, MNRAS, 371, 1983 
Martmez-Delgado D., Zinn R., Carrera R., Gallart C, 

2002, ApJ, 573, L19 
Martmez-Delgado D., Dinescu D.I., Zinn R., Tutsoff A., 
Cote P., Boyarchuck A., 2004, in Prada F., Martmez- 
Delgado D., Mahoney T.J., eds, ASP Conf. Ser. 327, Satel- 
lites and Tidal Streams. Astron. Soc. Pac, San Francisco, 
p. 255 

McConnachie A.W., Irwin M.J., Ferguson A.M.N. , Ibata 

R.A., Lewis G.F., Tanvir N., 2005, MNRAS, 356, 979 
McConnachie A.W., et al., 2009, Nature, 461, 66 
McConnachie A.W., Ferguson A.M.N. , Irwin M.J., Dubin- 

ski J., Widrow L.M., Dotter A., Ibata R., Lewis G.F., 

2010, ApJ, 723, 1038 
Olszewski E.W., Saha A., Knezek P., Subramaniam A., de 

Boer T., Seitzer P., 2009, AJ, 138, 1570 
Perina S., Bellazzini M., Buzzoni A., Cacciari C, Federici 

L., Fusi Pecci F., Galleti S., 2012, A&A, 546, 31 
Richardson J.C., et al., 2011, ApJ, 732, 76 
Schlegel D.J., Finkbeiner D.P., Davis M., 1998, ApJ, 500, 

525 

Searle L., Zinn R., 1978, ApJ, 225, 357 
Sohn Y.-J., et al., 2003, AJ, 126, 803 

Sollima A., Gratton R.G., Carballo-Bello J. A., Martmez- 
Delgado D., Carretta E., Bragaglia A., Lucatello S., 
Penarrubia J., 2012, MNRAS, in press (|arXiv:1207.6263P 

Stetson P.B., 1987, PASP, 99, 191 

Stetson P.B., et al., 1999, AJ, 117, 247 

Stoughton C, et al., 2002, AJ, 123, 485 

Tonry J., Davis M., 1979, AJ, 84, 1511 

Zinn R., 1993, in Smith G.H., Brodie J.P., eds, ASP Conf. 
Ser. 48, The Globular Cluster-Galaxy Connection. As- 
tron. Soc. Pac, San Francisco, p. 38 



© 2012 RAS, MNRAS 000, [THU 



